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We describe the singularities in the averaged density of 
states and the corresponding statistics of the energy levels in 
two- (2D) and three-dimensional (3D) chiral symmetric and 
time-reversal invariant disordered systems, realized in bipar- 
tite lattices with real off-diagonal disorder. For off-diagonal 
disorder of zero mean we obtain a singular density of states 
in 2D which becomes much less pronounced in 3D, while the 
level-statistics can be described by semi-Poisson distribution 
with mostly critical fractal states in 2D and Wigner surmise 
with mostly delocalized states in 3D. For logarithmic off- 
diagonal disorder of large strength we find indistinguishable 
behavior from ordinary disorder with strong localization in 
any dimension but in addition one-dimensional 1/\E\ Dyson- 
like asymptotic spectral singularities. The off-diagonal dis- 
order is also shown to enhance the propagation of two inter- 
acting particles similarly to systems with diagonal disorder. 
Although disordered models with chiral symmetry differ from 
non-chiral ones due to the presence of spectral singularities, 
both share the same qualitative localization properties except 
at the chiral symmetry point E = which is critical. 

Pacs numbers: 71.23.An; 73.20.Jc, 74.40.+k 



I. INTRODUCTION 



In the last couple of decades the effect of various sym- 
metries in the scaling theory Q of Anderson localiza- 
tion has been explored, both theoretically and experi- 
mentally, leading to fundamental results For exam- 
ple, in two-dimensional (2D) disordered systems breaking 
of the symmetry under time-reversal leads to the Quan- 
tum Hall Effect where the presence of magnetic field is 
responsible for critical delocalized states at the center 
of each Landau band which carry the Hall current ||, 
while breaking of symmetry under spin-rotation due to 
spin-orbit coupling is believed to lead to a metallic phase 
even in 2D Q. Recently, strong interest has focused on 
other basic symmetries, such as the chiral or particle- 
hole symmetry. In addition to the three standard uni- 
versality classes (orthogonal, unitary and symplectic) the 
presence of chiral symmetry gives three extra chiral en- 
sembles while four similar but distinct ensembles 
are realized for particle-hole symmetry in normal metal- 
superconducting systems, increasing their total number 
to ten The prerequisite for chiral symmetry is the 
presence of bipartite lattice structure with two intercon- 
nected sublattices and the simplest case where it can be 



realized in disordered systems is real random hopping 
between nearest- neighbours ||. This is also known as 
off-diagonal disorder and it is worth studying the spec- 
tral and localization properties of such systems in order 
to see the effect of chirality in 2D, 3D. The presence of 
chiral symmetry in disordered systems somehow simpli- 
fies both analytical and numerical calculations so that 
more is known about Anderson localization for the chiral 
symmerty classes in one, quasi-one and two dimensions il- 
lustrated for quasi-one dimension in || . In the rest of the 
paper we will rely on numerical tools in order to explore 
quantities that have not been computed analytically so 
far and also for quantities where the analytical predic- 
tions are competing, such as the 2D density of states. 

The chiral symmetry is present in simple nearest- 
neighbor models with off-diagonal disorder defined on bi- 
partite lattices which consist of two sublattices A and B, 
one connected to the other by random bonds B . This ge- 
ometric kind of symmetry, where the randomness strictly 
connects sites of two different sublattices, is different 
from other intrinsic symmetries of the Hamiltonian such 
as time-reversal and spin-rotation. It can also appear as 
an electron-hole transformation in half filled many body 
systems |7],|lO| . It is well-known that a consequence of the 
chiral symmetry is a distinct special point in the energy 
spectrum, for off-diagonal disorder in square and cubic 
lattices this is the band center E = with the eigenval- 
ues making up symmetric energy pairs {E, —E) around 
the zero mode E = 0. For an early observation of lo- 
calized chiral zero-modes see Ref . |ll| . For finite lattices 
the special chiral energy eigenvalue E = exists for any 
disorder configuration but only when the total number of 
sites is odd. Moreover, the corresponding E = wave- 
function is non-zero in one sublattice only and can be 
easily constructed p2 13 1. In fact, for any bipartite lat- 
tice with off-diagonal disorder only and ha (ub < tia) 
sites belonging to sublattice A (B) the number of lin- 
early independent E — states is exactly tia — hb , with 
their amplitudes being non-zero on the A sublattice and 
exactly zero on the B sublattice. The chiral E = state 
turns out to be neither localized nor extended but critical 
in all dimensions, also being very sensitive to the choice 
of boundary conditions Jl3] , reminiscent of critical states 
at the mobility edge of the Anderson transition . The 
main issue discussed in this paper is the singular behav- 
ior of the density of states in 2D, 3D, when approaching 
the special chiral symmetry energy E — where the lo- 
calization length also diverges Jq,ITq] . 
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The interest in this problem has been recently revived, 
due to the related problem of massless Dirac fermions in 
a random gauge fields in jlTj followed by |l8|,[l9|,|o| and 
also the Bogoliubov-de Gennes (BdG) Hamiltonians pl| ] 
which describe superconducting quasiparticles in mean- 
field theory. The Dirac fermions have very different zero 
disorder limit from the off-diagonal disorder case. In the 
tight-binding approximation the 2D fermions have a zero 
energy Fermi surface with many points forming a square 
while the Dirac fermions have a Fermi surface with four 
points only. In the former case the pure 2D density of 
states at E = has the well-known log-type van Hove 
singularity while in the Dirac case the density of states 
approaches zero at the band center . 

First, let us summarize what is known about the spec- 
tral and wavefunction anomalies of one-electron spec- 
tra in the presence of off-diagonal disorder in bipartite 
lattices described by real symmetric random Hamiltoni- 
ans. The ID random hopping problem in a bipartite 
lattice has a chiral spectrum with a strong l/\Ehi 3 \E\\ 
leading Dyson singularity of the density of states as E 
approaches zero, which is the first result obtained in 
disordered systems p2fl , even before the Anderson the- 
ory of localization was proposed p^] . The correspond- 
ing states of the energy spectrum are Anderson local- 
ized with their wavefunction amplitude decaying expo- 
nentially as IV'e^oMI ~ e _7r , where r is the distance 
and 7 the inverse localization length which depends on 
the disorder. This is not true at the chiral symmetry zero 
mode E = 0, where the corresponding wavefunction am- 
plitude which lies in one sublattice decays slower than 
exponential, as | - 0£; = o('")| ~ e _7v/? at distance r from 
its maximum where 7 also varies with disorder. This 
result for the typical E — wavefunction can be eas- 
ily derived since its log-amplitude In WpE=o(r)\ at site 
r executes a random-walk in space || leading to the 
exponential square root decay. In higher-dimensional, 
bipartite systems such as the square 2D or cubic 3D 
lattices the E = wavefunction was shown to display 
multifractal fluctuations with many scattered peaks and 
a disorder-dependent fractal dimension Moreover, 
an appropriate correlation function behaves logarithmi- 
cally which roughly implies an overall power-law decay 
1^=0(7")! ~ r~ v of the peak heights from the maximum 
peak and disorder dependent-??. It must be understood 
that such a power-law description is only a simplified pic- 
ture for an approximate decaying character of the chiral 
E = critical state. For chiral systems the density of 
states at the band center does not diverge only in ID. In 
the absence of time-reversal symmetry on a 2D bipartite 
lattice for weak disorder in the vicinity of the band center 

the divergence P2d{E) ~ \E\~~ 1 exp ^— c (in ^ ^ ^ was 

predicted by Gade J^] with a constant c and a universal 
exponent k = 1/2. If universality applies this singularity 
should also hold for the chiral-orthogonal ensemble. For 



strong disorder Mortunich, Damle and Huse |24| argued 
that the exponent k = 1/2 does not apply for a typi- 
cal density of states in the chiral orthogonal universality 
class and predict k = 2/3 instead. In |^5| the density of 
states for the 7r-flux phase is discussed via a weak disor- 
der approach. The analytical work of Q gives a wealth 
of exact results for the 7r-flux phase weakly perturbed 
by real nearest-neighbour hoppings. Recent work [ p7[ for 
random Dirac fermions instead of vanishing density of 
states shows an upturn close to E — 0, in agreement with 
the diverging dos of Gade ||. The multifractal analysis 
of |28| is compared to the analytical results of 20 1 . 



In this paper we shall investigate the consequences of 
the chiral symmetry for real symmetric random matrices 
in the orthogonal universality class by studying the spec- 
trum of 2D, 3D square and cubic tight binding lattices 
with off-diagonal disorder. In the studied systems apart 
from chirality the symmetries of time-reversal and spin- 
rotation are preserved as well. We shall determine the 
midband singularities in the density of states p{E) and 
from the nearest-level-statistics the localization proper- 
ties of the corresponding eigenstates. Our results demon- 
strate the effect of chirality for the unique off-diagonal 
disorder of mean zero: In 2D the nearest level-spacing 
distribution function P{S) turns out to be intermediate 
between Wigner and Poisson roughly described by the so- 
called semi-Poisson curve |l5| (for small-spacing S shows 
Wigner-like linear increase oc S and for large-spacing 
S shows Poisson-like exponential decrease oc exp{— S)), 
which indicates the presence of critical states in finite 2D 
disordered systems. In 3D our results for the eigenvalue 
spacing fluctuations with zero mean off-diagonal disorder 
are very close to the Wigner surmise which corresponds 
to extended states. For logarithmic off-diagonal disor- 
der which permits its strength to vary we see for very 
strong values that the distributions both in 2D, 3D be- 
come Poisson-like independently of the dimension and 
the states localize similarly to states of ordinary dis- 
ordered systems with broken chiral symmetry. Apart 
from clarifying the effect of chiral symmetry for non- 
interacting electron in disordered media our aim is also to 
extend our calculations to more than one electron by in- 
cluding Hubbard electron-electron interactions together 
with off-diagonal disorder, via the simplest possible two- 
electron Hamiltonian. This is done in order to discuss 
a recent conjecture concerning the decisive role of chi- 
ral symmetry for quantum transport properties in the 
presence of both disorder and interactions at half-filling, 
where enhancement of the conductivity with increasing 
disorder was recently obtained only for non-chiral sys- 
tems at half-filling |||. 

We present a numerical study of the spectrum and its 
fluctuations in squared and cubic bipartite lattices with 
disorder in the nearest-neighbor hoppings. We summa- 
rize our main questions: (i) What is the nature of the 
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spectral singularities in the presence of chiral symmetry? 
(ii) What is the corresponding level-statistics and the 
localization behavior for disordered systems when chi- 
ral symmetry is preserved? (iii) From the established 
spectral and localization behavior for systems with chi- 
ral or broken chiral symmetry can we conclude whether 
all states remain localized in 2D? (iv) How is the electron- 
electron interaction affected by off-diagonal disorder for 
two particles only? Apart from the theoretical interest 
in answering the above questions we emphasize that our 
study of localization with off-diagonal disorder can have 
a lot of applications. In particular, the presence or not 
of the chiral symmetry could be useful for understanding 
various properties of several realistic 2D systems, such 
as current-currying states in the quantum Hall effect 
||, quasiparticles in dirty superconductors 27 1, states 
in semiconductor quantum wells, high mobility silicon 
MOSFETs M, etc. 



II. THE DISORDERED HAMILTONIAN WITH 
CHIRAL SYMMETRY 



We consider a tight-binding model Hamiltonian de- 
fined in bipartite lattices with random nearest-neighbor 
hoppings of the off-diagonal block symmetric structure 



H = ^(^ C J Cj - + H.c. 







Ha,b 




(1) 



where the sum is taken over all bonds (ij) where i,j 
denote nearest neighbor lattice sites which belong to dif- 
ferent sublattices, Cj is the annihilation operator of an 
electron on site i and Uj the nearest-neighbor hopping 
integrals which are real independent random variables 
satisfying a box probability distribution with mean zero 
and arbitrary variance. This is equivalent to choosing 
random hoppings Uj from a box distribution P{tij) = 
— , for — ^ < ta < It of zero mean and width w. The 

w ' 2 — L J — 2 

scaling U : j/w is sufficient to make the disorder indepen- 
dent of u;, due to absence of other energy scale in the 
Hamiltonian. Thus, by increasing w the energy band be- 
comes wider since it amounts to a rigid rescaling of the 
eigenvalues only (the eigenstates remain the same for any 
w) and this is a unique form of zero mean off-diagonal 
disorder which allows certain comparison with analyti- 
cal results to be made. The alternative choice is disor- 
der of mean zero and width W but for the ln^.j. This 
logarithmic disorder distribution is PF-dependent, can be 
broad becoming arbitrarily strong and guarantees always 
positive hoppings U,j- It is commonly used in ID where 
the real hoppings t^.j chosen from a logarithmic range 
must be always positive and can be widely varied with 
typical W, which is believed to be a good measure of off- 
diagonal disorder more generally. We emphasize that in 



our study the random distribution of the hoppings Uj is 
the only disorder ingredient in the Hamiltonian H , since 
we completely ignore diagonal disorder by setting all site 
energies equal to zero. The diagonal disorder, somehow, 
connects each site to itself breaking the chiral symme- 
try. In our studies we consider finite L x L square and 
Lx Lx L cubic lattices with various choices of boundary 
conditions, taking care to preserve the bipartite structure 
and chirality of H . 

In Eq. (I) the matrix Ha,b contains matrix elements 
which connect the A, B interconnected sublattices defin- 
ing the two A, B bases and due to chiral symmetry the 
Hamiltonian transforms as 



H = -cr 3 H<7 3 , cr 3 



1 

—I 



(2) 



which means anticommuting H and 03 . One notices that 
03 behaves like the 75 familiar from applications of gauge 
theories fL7fl . Moreover, the presence of chiral symmetry 
allows to reduce the size of the corresponding matrices 
to half by diagonalizing H 2 instead of H obtaining the 
eigenvalues squared. Due to the interconnected by the 
random hoppings sublattice structure of the studied lat- 
tices Jl2| the eigenvalues appear in pairs E, —E around 
E = so that the spectral density is strictly an even 
function p(E) = p{—E) of the energy E. The eigen- 
value pairs E,—E around E = have simply related 
eigenstates ipsir) = J2 n ,i,... Vv.-K ^ ■■■) and ^P-E^r) = 
J2n,i,...(^ 1 ) n+l+ '"^i,i,...\n,l,...) since Hip = Eip from 
Eq. (2) gives Ha 3 ij) = —Ecf 3 ip so that the unitary trans- 
formation CT3 changes the wavefunction amplitude of the 
—E state on every other site, say on the sites of sub- 
lattice B, keeping the amplitudes of the —E state the 
same as that of E on sublattice A. Moreover, the E = 
state when it exists (for odd total number of sites) has 
amplitude only on A sublattice and is strictly zero on B 
sublattice. 

In order to measure the anomalies in the spectral den- 
sity and their corresponding fluctuations we obtained 
the eigenvalues of finite Hamiltonian matrices of size 
L d x L d , d = 2,3 for 2D, 3D, respectively. The very 
sparse matrices are obtained from Eq. (1) by distribut- 
ing either t% j uniquely or lnijj, with various strengths 
W to cover a wider range of off-diagonal disorder. We 
have used standard diagonalization algorithms, includ- 
ing Gaussian elimination for computing the integrated 
density of states (idos) by employing eigenvalue count- 
ing theorems at fixed E j3l[] and the Lanczos algorithm 
for computing eigenvalues within an energy range |32j. 
Moreover, we have considered all possible boundary con- 
ditions |lq] , such as periodic in all directions, periodic in 
one and hard wall in the other or hard wall in all direc- 
tions, always preserving the chiral symmetry. It must be 
also noted that crucial even-odd size effects Bp3| may 
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arise in this problem, for example, with periodic bound- 
ary conditions the sites of the A sublattice are connected 
to the sites of the B sublattice preserving chirality only 
for even linear size L. The periodic boundary conditions 
join sites which no longer belong to different sublattices 
for odd L and the chiral symmetry is broken. For hard 
wall boundary conditions no such even-odd effects arise. 
We should also remember that the zero mode E = for 
chiral systems exists for odd sizes only. More details of 
the numerical methods can be found in Rcfs. pupa. 



III. ZERO MEAN OFF-DIAGONAL DISORDER 
A. The singularities in the dos 

In one dimension (ID) we cannot have zero mean off- 
diagonal disorder due to the finite probability of finding 
a zero bond which breaks the chain. Outside one dimen- 
sion for this kind of disorder some rigorous results are 
known, such as the Gade || and Motrunich et al p^ j 
singularities. Earlier approximate results were derived 
with a jt - expansion |m| since the studied model is an 
example of the ensemble introduced in Ref. Jl(| . The ex- 
trapolation of the one-loop results near E = |32| gives 
a logarithmic singularity in 2D dos 



p 2D {E) = 0.299577. 



2tt 2 ln \E\ 



0{\E\), (3) 



similar to that of the pure system where the constant 
term 0.299577 is replaced by ^(lnlG + (1/2) ln2) « 
0.1580186 |]32| , |3~4j| . In the two-loop approximation the sit- 
uation is less clear, for example, the corresponding term 
is divergent, proportional to (ln pjy) 2 j which might imply 
a power-law p 2 d{E) cx |E| _( ^ spectral singularity. The 
small \E\ expression obtained in 2D by Gade Q com- 
bines the strong 1/\E\ divergence (Dyson- type) with a 
slowly decaying exponential of the logarithm of E to the 
power k, written as their product 



p 2D (E) = ^- exp (-c 2 (ln-L 



k = 1/2, (4) 



where ci,C2 are model-dependent constants. The im- 
proved strong disorder expression obtained by Motrunich 
etal |Q involves a different exponent for the logarithm 
K = 2/3. In 3D the one- loop result gives a square-root 
singularity which remains for higher loops implying finite 
dos at E — with 



dos p(E) close to E = 0. The random hoppings tjj are 
uniformly distributed for example within [—1/2, +1/2] 
since, of course, this specific range is irrelevant for local- 
ization. In order to determine the spectral singularities 
two possible fits are tried: 

(i) A power-law fit for the idos Af(E) cx \E\ 1 ^ 4 ' which 
implies for the dos the power-law divergence p(E) cx 

In Fig. 1(a) numerical data for the idos Af(E) 
which are obtained from eigenvalue-counting techniques 
allow to compute the exponent cj) for squared and cu- 
bic lattice with lattice sizes up to 600 x 600 in 2D and 
30 x 30 x 30 in 3D. In the log-log plots the data are seen 
to collapse to power-laws which define a small exponent 
<\> k, 0.25 in 2D and cj) w 0.07 in 3D. These weak power- 
law divergencies of the dos with off-diagonal disorder are 
certainly different from the expected log-type singulari- 
ties of non-random systems in pure 2D and the constant 
dos in pure 3D (see Fig. 1). 

(ii) We have also tried an alternative fit of the data in 
Fig. 1(b) by plotting directly p(E) vs. E in a semi- 
log plot. In 2D the data can be seen to be reason- 
ably explained by the analytical form of Gade H and 
Motrunich etal Q. The fit of Gade with n = 1/2 gave 
ci pa 21.5, C2 ~ 3.77 and the fit of Motrunich etal with 
k = 2/3 gave c\ ss 2.61, c 2 ~ 2.14. We have also tried an- 
other fit of the data which gave a slightly weaker power- 
law in the exponential with the value n « 0.91 and con- 
stants ci pa 0.47, C2 ~ 1.04. We also see that the one-loop 
result of Eq. (5) in 3D is clearly inappropriate to fit the 
small- IE 1 1 data. Other fits can be also tried such as the 
fit of Ref. H which instead of 1/|-E7j from Eq. (4) in- 
volves the square-root l/y^Ef leading divergence. The 
above fits (i), (ii) verify that the presence of zero mean 
off-diagonal disorder leads to a weakly singular dos at 
least in 2D and a much weaker singularity, or a simple 
increase of the dos close to E = 0, in 3D. 



p 3D (E) = 0.400494. 



0(\E\). (5) 




In Fig. 1 numerical results are presented for the aver- 
aged integrated dos (idos) M(E) — dE'p(E') and the 
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Fig. 1. (a) Log-log plot of the 2D and 3D integrated density 
of states Af[E), which measures the number of eigenvalues from 
to E, with power-law fits of the numerical data for zero mean off- 
diagonal disorder. The broken lines give the corresponding idos for 
the pure systems while the inset shows linear plots of the densities 
in the random case which display singularities at E = 0. (b) The 
2D, 3D density of states p(E) with the corresponding analytical 
forms of one loop of Gade [B with k = 1/2, of Motrunich et 
al Ei with k = 2/3 and our fit with k « 0.91. 



B. The level-statistics 



In order to check whether chirality has some effect on 
the localization properties we proceed with the compu- 
tation of spectral fluctuations via the nearest-level P(S) 
distribution function. First we "unfold" the spectrum 
by making local rescaling of the energy scales to make 
the mean level spacing A(E) oc l/p(E) constant, equal 
to one, throughout the spectrum. This is particularly 
important for a singular dos since the mean level spac- 
ing tends to zero at a singular point. Instead of the raw 
levels Ei after "unfolding" the spectrum consists of the 
"unfolded" levels 



£i=Af av (Ei), i = I, 2, ...,7V, 



(6) 



obtained from the averaged integrated spectral density 
Afav(E). For the "unfolded" levels £i the average nearest- 
level spacing becomes equal to one. Moreover, for the 
localization properties of chiral systems we should dis- 
tinguish between the majority of levels in the spectrum 
E 7^ and the state at the chiral point E = 0. 

In Fig. 2 we present the obtained 2D nearest-level dis- 
tribution function P(S) for zero mean off-diagonal disor- 
der, various sizes L and boundary conditions (BC), us- 
ing the majority of the energy levels in the band except 
those close to the band center E — and the tails. In 
2D the results for various L seem to lie in-between the 



Wigner surmise (extended states) and the Poisson distri- 
bution (localized states), described by the intermediate 
semi- Poisson distribution P(S) — 4S*exp(— 25) p5l. For 
small spacing P(S) is Wigner-like and for large spacing 
the tails seen in the figure insets are Poisson-like. We 
have also performed averages of P(S) over BC which can 
be also seen to be roughly described by the semi-Poisson 
curve Jl5| . It must be stressed that this distribution ap- 
pears only for intermediate sizes with L smaller or close 
to the localization length £ |T^,|T|,^5| in the considered 
energy window. The obtained results also agree with Ref. 
p6[ where it is shown that the majority of the states for 
finite 2D disordered systems are multifractal for certain 
sizes L below £ (£<£). 
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"unfolding" over energy by using the well-known analyt- 
ical form N av {E) = (W 2 /9)/(ln\E\) 2 in ID @, work- 
ing instead with the "unfolded" levels Si oc 1 / (ln\Ei\) 2 , 
i = 1, 2, N. We observe that the spacing distribution 
of the first positive eigenvalue from the zeroth shown 
in (a) is distinguished from the other spacing distribu- 
tions (b), (c) and clearly shows a higher degree of level- 
repulsion. In this case the localization length £ can be 
also obtained for every energy E from the expression 
£(E) = (12/W 2 ) \ln\E\\ (e.g. see @). In Fig. 3 the 
considered size L is always larger than £ but P(S) which 
depends only on the ratio L/£ reaches the Poisson dis- 
tribution only asymptotically, for L much larger than £. 
We emphasise again that for E close to the chiral en- 
ergy E = (where £ is larger) the Poisson law for ID 
localization is seen to be reached more slowly. 



P(S) 




Fig. 2(a),(b),(c),(d). The 2D nearest level-spacing distri- 
bution function P(S) for zero mean off-diagonal disorder, various 
sizes L and three possible boundary conditions, such as periodic 
(PBC), hard wall (HWBC), periodic+hard wall and also aver- 
age over them. The dashed line is the Wigner surmise P(S) = 
(7r/2)S'exp(-(7r/4)5 2 ), the dotted line is the Poisson law P(S) = 
exp(— S) and the solid line is the intermediate critical semi-Poisson 
distribution P(S) = 4Sexp(— 25"). In the insets we display log-log 
plots for the corresponding tails of P(S). For L = 20, 30, 40, 50, 70 
the number of runs is 5000, 6000, 10000, 3555, the eigenvalues were 
located in the energy band between [2, 6], for w = 10 where the full 
energy range is [0, 13], and their total number was a few million. 

In order to demonstrate the effect of stronger level- 
repulsion close to the chiral E — state we choose as 
an example a ID chiral system with off-diagonal disor- 
der for L > £ and odd number of sites for the E = 
state to exist. In Fig. 3 we plot the corresponding 
ID spacing distribution P(S) of the spacings S defined 
from the differences between the first from the zeroth 
(Ei — 0), the second from the first (E 2 — £q.) and the 
eleventh from the tenth (En — E w ) eigenvalues for log- 
arithmic disorder of W = 1. This involves "unfolding" 
over the random statistical ensemble since the statistics 
of each spacing S is studied for many random configu- 
rations at fixed energy. We have additionally performed 



Fig. 3. The ID nearest level-spacing distribution function P(S) 
of (a) the first from the zeroth (with localization length £ 105), 
(b) the second from the first (£ 75) and (c) the eleventh from 
the tenth (£ «3 36) eigenvalues for logarithmic off-diagonal disorder 
W = 1 and larger size L = 1001. The continuous line is the 
asymptotic Poisson law P(S) = exp(— S) valid for localized states 
in the limit L>$ and the total number of runs was 500000. Note 
the double "unfolding" performed over both disorder and energy 
since we are very close to the singular dos. 

In Fig. 4(a), (b) we return to the 2D zero mean case and 
plot the probability amplitude distribution of two typi- 
cal multifractal wavefunctions at energies E ^ 0. They 
are seen to spread over space with fractal characteristics 
described by the fractal dimension D 2 and display many 
scattered peaks (scars) . It must be pointed out for £ ^ 
the states are multifractal only for L < £ where £ is their 
localization length, while the states at E = where £ 
diverges are described by multifractal dimensions valid 
for any size L ]3q] . In order to determine £ instead of 
the exponential decay of the wavefunction amplitude on 
very long strips, alternatively we can count the number 
of sites where significant amplitude exists via the inverse 
participation ratio (IPR) ]i3| , p5[ . The computed IPR is 
a rough measure of £ since it favors only the large ampli- 
tudes rather than the exponential tails of the amplitude 
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distributions. The results for the average IPR in 2D are 
displayed vs. energy E in Fig. 4(c), (d) for various sizes 
L with hard wall and periodic BC, respectively. They al- 
low to estimate that localization is not yet reached for the 
studied energy region with the considered sizes L lower 
or equal to the localization length £. For example, in 
Fig. 4(c), (d) we can see that for the largest size L = 100 
the mean IPR is close to converge, which implies that £ 
is roughly of the same order of magnitude [L sa £). In 
Fig. 4(c), (d) localization should be reached only when 
the curves of the mean IPR for two different sizes coin- 
cide. We can also speculate that the semi-Poisson ap- 
pears because in the considered energy range for P(S) in 
Fig. 2 we encounter a crossover situation where £ which 
does not vary a lot with E obeys L < £, that is L roughly 
smaller or equal to £. The intermediate semi-Poisson dis- 
tribution found in 2D apart from intimately connected to 
the underlying multifractality of the wavefunctions (for 
L < £) can be regarded as a crossover distribution with 
L w £. We emphasize that the adopted energy range lies 
in the middle between E — and the tail with the chosen 
sizes L either smaller or close to the wavefunction radius 
estimated by the corresponding £. 
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Fig. 4(a), (b). The amplitude distributions of two typical 
eigenstates for zero mean off-diagonal disorder in 2D away from 
the band center JJ/0. The states are multifractal only on scales 
below the localization L < £ length . The behavior of the E = 
state when it exists is different since this state is critical for any 
size with fractal dimension D2 ranging from the space dimension 
D2 = d for weak off-diagonal disorder to D2 = for strong off- 
diagonal disorder Jlijj . (c),(d). The corresponding average inverse 
participation ratio {IPR) vs. energy E for various chiral even sys- 
tems with hard wall and periodic BC, respectively. Note the values 
of (IPR) are lower for PBC where the states are slightly more ex- 
tended in comparison. For finite systems from IPR <x L~ D-2 near 
the minimum of the curves we find D2 ~ 0.8, slightly higher D2 
near E = and much lower D2 (localized states) for large E. The 
absence of saturation of the curves implies that the localization 
length is larger than the system for the adopted sizes (L < £). In 
the bacground with grey squares we see the actual values of IPR 
for the size L = 100 only. They display enormous fluctuations and 
are used for computing the corresponding average for L = 100. 

The obtained semi-Poisson P(S) distribution and the 
multifractality for finite size in 2D concerns states only 
within the band (E ^ 0) for finite L < £. For the zero 
mean off-diagonal disorder the chosen sizes L are smaller 
or equal to the localization length £. However, for larger 
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L ^> £ (we cannot reach such L in our 2D calculations) 
the data from Fig. 2 should eventually approach the Pois- 
son curve which indicates only localized states in 2D. The 
very slow approach to the localized Poisson limit for in- 
finite size L — ► oo is expected here except at the band 
center E = where £ also diverges [ fl6||35| . In 3D the 
obtained P(S) for this zero mean off-diagonal disorder is 
shown in Fig. 5 to lie very close to the Wigner surmise, 
a result compatible with the existence of mostly delocal- 
ized states in 3D. This implies the asymptotic approach 
towards Wigner for this kind of off-diagonal disorder. In 
other words, the obtained behavior exhibited by P(S) for 
E suggests critical behavior only for finite systems 
in 2D and delocalized diffusive behavior for any size in 
3D. Therefore, the localization properties of chiral sys- 
tems outside the chiral energy E — in any dimension is 
similar to that of ordinary disordered systems with the 
corresponding localization length £. 




Fig. 5. The 3D nearest level-spacing distribution function 
P{S) for the zero mean off-diagonal disorder, various system 
sizes L with the dashed line giving the Wigner surmise P(S) = 
(7r/2)Sexp(— (7r/4)5 2 ). The sizes L = 5,10,15 shown are diago- 
nalized for 5000, 6000, 2000 configurations and a few million total 
number of eigenvalues was considered in each case. 

We have also checked the probability distributions for 
the first few eigenvalues close to E — by employing 
many random configurations (realizations of disorder) as 
done in [|35f . This demonstrates nicely the level-repulsion 
(not shown) by using data for lattices of even number of 
sites with periodic BC and odd- number of sites with hard 
wall BC, which preserve chiral symmetry. Although chi- 
rality is present in both cases the E = state exists 
only for the later odd case. We have observed the in- 
crease of the eigenvalue density p(E) close to E = for 
both even and odd sizes but could not draw any con- 
clusions for the statistical behavior of the energy levels 
since energy level-repulsion can be examined only after 
"unfolding" . For the even size L systems the two eigen- 



values closest to E = (the lowest positive one and the 
one with exactly the same energy but minus sign) al- 
ways show level-repulsion due to chiral symmetry. The 
level- repulsion between the two energies E\, —E\ for the 
smallest possible E\ is expected since these energies are 
always separated by the chiral symmetry point E = 0. 
On the other hand, the effect of level-repulsion should 
also be enhanced for states close to E = where the 
localization length £ diverges (longer £ implies stronger 
level- repulsion). For odd size L with chiral symmetry the 
E = eigenvalue exists for every random configuration. 
The conclusion which can be drawn for this chapter is 
that scaling of level-statistics for E ^ chiral systems is 
determined by the ratio of the system size L over the cor- 
responding localization length £, similarly to what hap- 
pens for ordinary disordered systems and a critical state 
exists at E = 0. 



IV. LOGARITHMIC OFF-DIAGONAL DISORDER 



A. The 1/\E\ leading singularity 



In 2D the states away from the band center [E ^ 0) 
are also weakly multifractal on scales below the localiza- 
tion length. The discontinuous behavior claimed for the 
fractal dimension D2 at the band center in is proba- 
bly due to the different boundary conditions used for the 
E = and E ^ states. At E = where the local- 
ization length diverges the corresponding state is multi- 
fractal for every length scale with a fractal dimension D2 
dependent on the strength of the logarithmic off-diagonal 
disorder W , ranging from 2 to for very weak and very 
strong off-diagonal disorder, respectively. For strong off- 
diagonal disorder of width W for the In tij the observed 
accumulation of levels close to E = behaves differently. 
The data shown in Fig. 6 for 2D, 3D can be fitted to 
power-law singularities for the idos 



Af(E) cx l^l 1 ^ <*p(E) oc \E\ 



(7) 



with the exponent equal to <f>. The AT(E) and the fitted 
exponents <f> are shown in Fig. 6 where the 1/\E\ Dyson 
singularity which implies (j> = 1 is approached when W in- 
creases. This indicates that the behavior depends on the 
strength W so that for very strongly disordered medium 
the density of states should approach its ID Dyson limit, 
as predicted in J3lj. This is also in agreement with the 
forms of Gade and Motrunich et al forms having very low 
constant C2 in Eq. (4). 
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Fig. 6. The exponent <j> for the power-law fit of the density 
of state p(E) oc \E\~ ' is shown in 2D, 3D to approach the ID 
Dyson singularity limit (power-law with cf> = 1). In the inset the 
corresponding idos for various values of the strength of logaritmic 
off-diagonal disorder is shown. Eigenvalue counting algorithms are 
used for sizes L = 200, 400 and L = 30 in 2D, 3D, respectively. 



B. The tridiagonalization method 

The defined chiral 2D or 3D Hamiltonian H can be 
mapped onto a one-dimensional semi-infinite chain via 
the so-called tridiagonalization scheme • The spec- 

tral singularities and the degree of localization is then 
reflected on the statistical properties of the created tridi- 
agonal matrix. This is a version of the Lanczos algorithm 
p8| which is normally a procedure for finding eigenvalues 
and eigenvectors of large sparse matrices with most ma- 
trix elements zero. The Lanczos method like any other 
diagonalization method reduces the matrix into a tridi- 
agonal form proceeding as follows: one starts with a nor- 
malized starting state $i (e.g. corresponding to a central 
lattice site) and by operating successively with H gen- 
erates the orthonormalized states $„, via the recursion 
formula 



/3 n $„ = {H- a„_i)$„_i - ^„_i$„_ 2 



(8) 



n = 2, ...,7V and (3\ = 0. The new states $ n ,n = 
2, N are linear combinations of the original lattice ba- 
sis set, e.g. in 2D this set consists of the normalized 
states |ni,n2), n\, n 2 = 0, ±1, ±2, ±3, .... In the new 
$i, $2, ■•■ basis H is represented by a semi-infinite tridi- 
agonal matrix where the diagonal matrix elements are 
a n , n = 1, 2, N and the off-diagonal /3„, n = 2, 3, N. 
The tridiagonalization preserves chiral symmetry so that 
a n = 0,n = 1,2, N for chiral systems. In that case if 
we start from the central site corresponding to the nor- 
malized state $i = 1 0,0) from Eq. (8) with zero a n 
we have $ 2 = H$i//3 2 , $3 = H^ 2 //3s, - (/3 2 /%)$i, 



etc. and the (3 n ,n = 2,3,... are obtained at each step 
by normalizing the states $2, $3, Since H is a near- 
est neighbor Hamiltonian acting on |ni,ri2) we can ar- 
range the lattice into "shells" . The state & n +i includes 
all the lattice sites up to nth "shell" (for zero disor- 
der only the sites of the nth "shell") which in 2D con- 
sists of the lattice sites m, n 2 = 0, ±1, ±2, ±3, ... so that 
n = \ni\ + |n 2 |, starting from $1 = |0,0). As an exam- 
ple, for the periodic 2D Hamiltonian one computes f3 2 — 
2,/3 3 = \/5, 04 = ^76/20, etc, approaching /3„— 00 = 2. 
The basis states for the corresponding tridiagonal matrix 
are$! = |0,0), $ 2 = (|1, 0) + 10, 1) + | -1, 0) + |0, -l))//3 2 , 
$3 = (|2,0) + |0,2) + |-2,0) + |0,-2)+2|l,l)+2|l,-l) + 
2|-l,l> + 2|-l,-l»/(ft03), etc. 

The advantage of the tridiagonalization method is 
economy of storage since we can save only the a's and 
/3's as we go along. The weakness of the method is loss of 
orthogonality due to rounding errors (although for sparse 
matrices the damage is smaller). In 2D or 3D if we start 
from a central site the matrix consists of blocks around it 
so that each block corresponds to a "shell" . A site which 
can be reached from the central site (0th "shell") in a 
minimum of n steps from the origin is in the "shell" n. 
The 2D or 3D matrix is block-tridiagonal with blocks cor- 
responding to "shells" starting from the 0th "shell" with 
$1 and we can work with N iterations, exactly equal 
to the number of "shells", to determine $2, $3, ■■■,&n+i 
and f3 2 , /?3, Pn+i for the 1,2, ...,iVth "shells" around 
the central site. From the obtained ID semi-infinite tridi- 
agonal matrix one can conveniently compute the local 
density of states on the central site $1 which has its first 
27V moments exactly equal to the spectral density of the 
original 2D or 3D matrix. It must be stressed that the 
use of the Lanczos algorithm nowadays is different than 
that of the tridiagonalization method. Its wide use is 
directed, instead, towards the computation of specific 
eigenvalues of large sparse matrices with the requested 
accuracy. This can be achieved by increasing the num- 
ber of iterations to much larger values than the number 
of "shells" or "blocks" considered in tridiagonalization. 

For zero disorder or disorder which preserves chiral 
symmetry connecting A and B sublattices a 2D or 3D 
tight binding system maps to a chain with all a n = 0. 
For chiral disorder the mapping onto a ID chain in- 
volves only off-diagonal disorder but inhomogenous type 
(n— dependent) I^Jjlf . It is only for very strong off- 
diagonal disorder where the disordered chain resulting 
from the tridiagonalization becomes homogenously off- 
diagonal disordered. On the basis of our numerical re- 
sults we argue that although this conjecture holds the 
suggested inhomogenous decay in the case of not so 
strong disorder is absent as well (or it is not of power-law 
type). 

Since for chiral systems the ID semi-infinite chain is 
also chiral with off-diagonal disorder only the well-known 
ID results could map to higher dimensions. The resulting 
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ID chain from the tridiagonalization scheme has inho- 
mogenous off-diagonal disorder with hopping matrix ele- 
ments (3 n chosen from an site n-dependent random distri- 
bution of mean {(3 n ) and variance var(/?„) = ((3%) — (Pn) 2 - 
Under the approximation that the states are ex- 

pressed as linear combinations of |ni,n 2 ) belonging to 
the nth "shell" only, which is valid for pure systems and 
claimed to hold for weak off-diagonal disorder, in |37]] 
a power-law decay of the variance vs. n was predicted. 
This gave ip(r) oc e -1 ^^ decay of the E = chiral 
state in every dimension higher than one. Our numerics 
show that some of the arguments which led to this re- 
sult are not valid since for weak off-diagonal disorder the 
variance vs. n does not decay as a power-law, but one 
eventually obtains indications of saturation for large n 
(although the data for smaller n could be fitted to power 
laws). For strong off-diagonal disorder the obtained sat- 
uration implies a constant var(/3„) independent of n and 
only in this case the results for ID off-diagonal disorder, 
such as the l/|i?| Dyson- type spectral singularities, carry 
through to higher dimensions. 
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Fig. 7(a). The parameters of the tridiagonal chain with mean 
(fin) and variance var(/3„) as a function of the chain index n for 



zero mean off-diagonal disorder and 5000 random configurations. 
Although decay of the variance is seen for small n, implying inho- 
mogenous disorder for the chain, the data point towards saturation 
at large n. (b).The parameters of the tridiagonal chain (/3„) and 
their variance var(/3„) as a function of the tridiagonal chain index 
n for various values W of logarithmic off-diagonal disorder. The 
saturation of the variance var(/3 n ) suggests a mapping of a higher- 
dimensional strongly off-diagonal disordered systems to a random 
chain. 

In order to check these predictions we are able to ex- 
tend the previous numerical calculations to higher num- 
ber of "shells" n starting from a central site <E>i (0th 
"shell"). For example, for the zero mean disorder in 2D 
we can go up to N = 1000 "shells" starting from the 
site at the origin and the total number of sites from the 
starting central site (0th "shell" ) up to the TVth "shell" is 
2N 2 + 2N + 1 . For the largest N used we always find sat- 
uration of the variance vs. n, even for the zero mean off- 
diagonal disorder (Fig. 7(a)) and certainly for the loga- 
rithmically strong off-diagonal disorder (Fig. 7(b)). The 
predicted inhomogenous power-law decay of p?]], which 
implies n-dependent variance, seems not to hold but such 
decay might exist only for finite sizes (small n). For large 
enough n for both zero mean and strong logarithmic off- 
diagonal disorder the mapping to the ID chain involves 
n-independent (homogenous) disorder. The above dis- 
cussion is not in favor of the square root of log expo- 
nential law predicted in the tridiagonalization of |57| and 
rather points towards the presence of 1/|-E| forms for the 
dos in higher d predicted in Gade and Motrunich et al. 
For strong off-diagonal disorder the n-independent vari- 
ance prediction with associated ID behavior in 2D, 3D 
holds more clearly and the 1/\E\ Dyson divergence should 
somehow appear. The data for the hoppings ((3 n ) and the 
variance var(/3„) vs. n are shown in Fig. 7(a) for zero 
mean and (b) strong logarithmic off-diagonal disorder. 
Although the saturation of (j3 n ) can be seen in all cases 
the var(/3„) seems to decay for the zero mean off-diagonal 
disorder. However, this is true only for small n and the 
var(/3„) is eventually seen having a tendency to saturate 
for larger n (Fig. 7(a)). For large W the saturation al- 
ready occurs for smaller n (Fig. 7(b)). In 3D we have 
obtained similar results. 



V. DYNAMICS OF TWO-INTERACTING 
PARTICLES WITH OFF-DIAGONAL DISORDER 



In Fig. 8(a) we observe the dynamics of one parti- 
cle for logarithmic off-diagonal disorder in ID. We find 
that the effect of off-diagonal disorder turns out to be 
similar to that of diagonal disorder, both deviating from 
the pure chain ballistic law a 2 (t) — 2t 2 (dashed line in 
Fig. 8(a)) in proportion to the added disorder. Next we 
turn to the interacting case since the chiral symmetry 
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(rather its absence) was recently shown to be crucial for 
the possibility of a metallic phase in 2D. For half- filled 
interacting systems the Mott gap turns out to be insensi- 
tive to weak disorder when chiral symmetry is present, so 
that a metallic phase in 2D could become possible only 
when the chiral symmetry is broken p9[ . 
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Fig. 8(a). Log-log plot of the mean-square displacement u 2 (t) 
as a function of log(t) in ID (note the double log(log(t)) in the 
x-axis so that 1, 6 correspond to i = 10, t = 10 6 , respectively, 
and the single log a 2 in the y-axis, all for base ten logarithms) for 
logarithmic off-diagonal disorder with W=0.l, W=0.5 and VF=1. 
The dashed line is the ballistic limit 2t 2 . In the inset results for di- 
agonal disorder from [— W/2, +W/2] are shown, (b) A log-log plot 
of the mean-square displacement cr 2 (t) for two particles as a func- 
tion of log(t) for logarithmic off-diagonal disorder with W=l and 
interaction strength U = 0, 1. The enhancement of the wavepacket 
in the presence of interaction as well as the smoother evolution 
with an approximate ultraslow diffusive law can be seen also for 
off-diagonal disorder. 

In order to check the effect of repulsive interactions in 
ID chain we have computed the mean square displace- 
ment from the formula 



a 2 (i) = (l/2) \ x fn 1 ,n 2 (t)\ 2 {n 1 -n 1 +n 2 -n 2 ) 2 , (9) 



with ^m,n2 00 the time dependent wavefunction for two 
particles at sites ni, n 2 in the presence of off-diagonal 
disorder and interaction U ||^,|l|, both starting at the 
same initial zeroth site. The mean electrons positions 
are n % = Y, ni ,n 2 \^n u n 2 {t)\ 2 ni, i = 1,2. The result 
for <r 2 (t) is shown in Fig. 8(b) for off-diagonal disorder 
without interaction (U = 0) and interaction of strength 
[7 = 1. We verify that the two-particle propagation in- 
creases with interaction, as for ordinary diagonal disor- 
der, where propagation is known to be enhanced in the 
presence of disorder when U is switched on 
increasing disorder we find that propagation decreases as 
expected. In the presence of both disorder and inter- 
action propagation is also seen to become smooth and 
the data obey an ultraslow diffusion law a 2 (t) oc log" (t) 
with a slope from the linear fit equal to a w 3 (Fig. 8(b)) 
when U = 1. Unfortunately, we cannot discuss the half- 
filled chiral symmetric limit since our case involves only 
low filling of two particles only. For U ^ our system 
clearly violates the chiral symmetry since it is no longer 
particle-hole symmetric due to U which is added in the 
matrix diagonal. More generally a bipartite lattice struc- 
ture with chiral symmetry was discussed in the absence 
of disorder pl[ | for general filling and dimension with in- 
teraction strength U. Their results hold for our chiral 
systems too since for negative U the spectrum turns out 
to be opposite to that for positive U while the energy 
E corresponds to the energy — E of the other model, re- 
spectively. This implies that for certain initial conditions 
the quantum evolution is the same for both positive and 
negative U @. 



VI. DISCUSSION 

We have investigated numerically 2D and 3D models 
of non-interacting fermions and ID model of two inter- 
acting fermions, all cases with random nearest-neighbour 
hopping in bipartite systems. They belong to the chiral 
orthogonal universality class where disorder respects chi- 
ral symmetry, manifesting itself in spectral singularities 
and certain anomalous localization behavior at E = 0. 
It should be pointed out that the effect of chiral symme- 
try in disordered systems was suspected more than two 
decades ago in studies of off-diagonal disorder [pi and 
about a decade ago in similar studies of BdG Hamilto- 
nians which describe dirty superconductors p2| . It was 
believed that exceptions from the scaling theory of local- 
ization might arise in these systems, such as spectral sin- 
gularities and a kind of derealization at the band center 
in an otherwise localized spectrum. Moreover, depend- 
ing on each author's viewpoint the critical behavior of 
the undoubtedly critical multifractal state at E = was 
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named either derealization |19| or power-law localiza- 
tion |0| , although both meant weakening of localization 
close to E = where the localization length diverges. 
The E = critical state can be fully described by mul- 
tifractal exponents via scaling the amplitude moments 
and its fractal dimensions, such as D2, vary strongly 
with the strength of off-diagonal disorder ranging from 
space filling extended states D 2 — > d for weak disorder 
to point-like localized states D 2 — > for strong disorder. 
This transition in the multifractal spectrum is known as 
a "freezing transition" |2C| . Below "freezing" computed 
mean and typical values agree but not above "freezing" 
where the means are determined by the tails of the cor- 
responding distributions. In 3D for weak off-diagonal 
disorder where most states with E ^ are extended the 
critical E — state also behaves as an extended state 
with D 2 ~ d. This picture also agrees with early trans- 
mission studies for particle-hole symmetric BdG Hamil- 
tonians which showed that quasi-extended states might 
exist Q or not depending on the chosen energy E. 
The present work summarizes all effects of chiral sym- 
metry such as spectral singularities in 2D for zero mean 
off-diagonal disorder which become strong 1/E Dyson- 
type singularities for strong off-diagonal disorder, while 
the corresponding level-statistics varies around the crit- 
ical semi-Poisson curve for finite systems of size smaller 
or equal than the localization length (L < £) suggesting 
that the E ^ states from multifractal become fully lo- 
calized for infinite size except the chiral E = energy 
state which is multifractal for any size, disorder and di- 
mension |44|j45|] • For strong off-diagonal disorder we find 
that the behavior for the rest of E ^ states is indistin- 
guishable from that of ordinary 2D disordered systems 
with the appropriate localization length. 

We find that the chiral symmetry realized in bipartite 
lattices when the disorder connects one sublattice to the 
other is characterized by enhanced dos and weakening of 
localization with a critical state at E = 0, having differ- 
ent spectral and localization properties from the E ^ 
states. This answers questions (i) and (ii) set in the in- 
troduction. Concerning question (iii) 2D derealization 
in this context seems to require interacting electrons and 
possibly broken chiral symmetry p9| . About question 
(iv), the effect of chiral disorder on the electron-electron 
interaction for two particles seems to be no different than 
that of diagonal disorder. The main message of section V 
is that off-diagonal disorder has the same effect as diag- 
onal disorder for coherent two interacting particle prop- 
agation, with the presence of ultraslow diffusion. Our 
study also offers answers to questions about the eigen- 
states. Outside the middle of the band (E 7^ 0) criticality 
is expected to occur only for finite systems [[36| . Instead, 
the E = state is critical with fractal exponents close to 
space filling for small W (extended states) and point-like 
for large W (localized states) KG] . This is similar to what 



one expects for ordinary Anderson localized states with 
diagonal disorder. 

In summary, we have presented numerical results for 
the density of states with off-diagonal disorder in order 
to test various theories which describe the observed spec- 
tral singularities, from Dyson to Gade and Motrunich 
etal forms. These results complement many analytical 
methods developed especially for unitary chiral systems. 
For zero mean off-diagonal disorder our main finding is 
that the form of the spectral singularity in 2D is stronger 
than the log-type singularity of the pure system and can 
be fitted either with a weak power-law or by the interme- 
diate Gade or Motrunich etal singularities and also their 
formula but with a new exponent k. In 3D the dos is a 
little enhanced at E — but when fitted to a power-law 
the obtained divergence, if any, is very weak. The corre- 
sponding level-statistics for finite systems in 2D is close 
to the critical semi-Poisson curve and in 3D the Wigner 
curve, implying critical and extended states for the ma- 
jority of states in the band for 2D, 3D, respectively. For 
logarithmically strong off-diagonal disorder our results 
are indistinguishable from ordinary strongly disordered 
systems except at the midband critical state E — where 
the chiral symmetry is responsible for anomalous dereal- 
ization. Emphasis is given in 2D where exceptions from 
generic localization are more pronounced in the form of 
large localization lengths and multifractality of all states 
for sizes below their localization length. However, for 
large scales a chiral symmetric disordered system behaves 
like an ordinary disorder system with only two visible 
differences: (i) Spectral singularities at the chiral energy 
E = which become 1/\E\ Dyson-like for strong disor- 
der or low dimensionality and (ii) the presence of mul- 
tifractality at the critical E = symmetry point which 
becomes close to extended for zero mean disorder and 
localized for strong disorder. At other energies the chiral 
disorder has the same effect as ordinary disorder. 
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